Unconventional sign-reversing superconductivity in LaFeAsOi_ x F 
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We argue that the newly discovered superconductivity in a nearly magnetic, Fe-based layered 
compound is unconventional and mediated by antiferromagnetic spin fluctuations, though different 
from the usual superexchange and specific to this compound. This resulting state is an example of 
extended s-wave pairing with a sign reversal of the order parameter between different Fermi surface 
sheets. The main role of doping in this scenario is to lower the density of states and suppress the 
pair-breaking ferromagnetic fluctuations. 



The discovery of superconductivity with T c > 26 K 
[l[ in a compound that contains doped Fe 2+ square lat- 
tice sheets raises immediate questions about the nature 
of the superconducting state and the pairing mechanism. 
A number of highly unusual properties suggest, even at 
this early stage, that the superconductivity is not con- 
ventional. We argue that not only is it unconventional, 
but that doped LaFeAsO represents the first example of 
multigap superconductivity with a discontinuous change 
of the order parameter (OP) phase between bands, a 
state discussed previously {e.g., Refs.@,|j|), but not yet 
observed in nature. 
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We suggest that superconductivity here is mediated 
by spin fluctuations (SF), as many believe is the case 
in cuprates, heavy fermion materials, or ruthenates. SF 
can only induce a triplet superconducting OP, or a singlet 
one that changes sign over the Fermi surface (FS). The 
latter condition often, but not always, dictates strong 
angular anisotropy of the OP (cf. d-wave). In our sce- 
nario it is satisfied despite full angular isotropy, since the 
sign changes between the two sets of FSs. Our model 
is also principally different from the convetional s-wave 
superconductivity discovered in MgB2: there the pairing 
interaction is attractive, in our case it is repulsive (but 
pairing thanks to the sign reversal). Finally, similar to 
d- or p-wave superconductivity our OP has a nearest- 
neighbor structure in the real space, thus reducing the 
Coulomb repulsion within the pair. 

We begin by arguing against conventional supercon- 
ductivity. The pure compound, LaFeAsO, is on the verge 
of a magnetic instability: it has a very high magnetic 
susceptibility [H and is strongly renormalized compared 
to density functional (DFT) calculations § . This renor- 
malization is higher than in any known conventional su- 
perconductor, including MgCNi 3 , where superconductiv- 
ity is believed to be strongly depressed by spin fluctua- 
tions. The susceptibility in the pure compound is large 
and upon doping with F grows even larger and becomes 
Curie- Weiss like. This suggests nearness to a critical 
point in the pure compound and non-trivial competition 
between different spin fluctuations (SF). Very strong elec- 
tron phonon interactions would be required to overcome 



FIG. 1: color online(a) Calculated band structure at x — 0.1 
near the Fermi Level, (b) Calculated Fermi surface at 10% F 
doping. Note that the only 3D parts are the far ends of the 
electron cylinders around M. The fully three-dimensional sur- 
face present in the undoped compound is suppressed beneath 
Ef by increased electron count. 



the destructive effects of such SF. We have calculated 
ab initio the electron-phonon spectral function, a 2 F(i 
and coupling, A, for the stoichiometric compound 
Some moderate coupling exists, mostly to As modes, but 
the total A appears to be ~ 0.2, with w/ ofl ~ 250K (u)i og is 
the logarithmically averaged boson frequency of Eliash- 
berg theory), which can in no way explain T c > 26 K. 

The calculated DFT Fermi surfaces [8] for undoped 
LaFeAsO consist of two small electron cylinders around 
the tetragonal M point, and two hole cylinders, plus a 
heavy 3D hole pocket around T. To study doping ef- 
fects, we performed full-potential calculations using the 
WIEN package in the virtual crystal approximation and 
the PBE GGA functional. The lattice parameters we 
took from experiment [l[ and optimized the internal po- 
sitions 0. In Fig. [D we show the bands near the Fermi 
level for x = 0.1 and the corresponding Fermi surface. 
As expected, the 3D pocket fills with electron doping (at 
x = 0.04-0.05) and the fermiology radically simplifies, 
leaving a highly 2D electronic structure with two heavy 
hole cylinders and two lighter (and larger) electron cylin- 
ders. 

This fermiology imposes strong constraints on the su- 
perconductivity. In particular, with the exception of 
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FIG. 2: color online Fermi surface formation upon backfolding 
of the large BZ corresponding to a simple Fe square lattice, 
(a) Real space: the four small unit cells of the Fe sublattice 
only (dashed) with the larger solid diamond of actual two- 
Fe unit cell. Dark and light circles indicate superexchange 
(checkerboard) ordering. The inset shows the spin density 
wave corresponding to X point SF. (b) Reciprocal space: the 
black square is the unfolded BZ, the blue diamond is the 
downfolded zone, the blue ellipse is the electron pocket from 
the Y point downfolded onto the X point (which is M in the 
small BZ). 

phonons, it is hard to identify pairing interactions with 
a strong k z dependence. Thus, states with strong varia- 
tions of the OP along k z are unlikely. An angular varia- 
tion of the OP in the xy plane is possible, but would re- 
quire an unrealistically strong q-dependence of the pair- 
ing interaction on the scale of the small Fermi surface 
size, and would also be extremely sensitive to impurities. 

On the other hand, the small Fermi surfaces are readily 
compatible with a pairing state with weak variations of 
the OP within the sheets, but a ir phase shift between 
electron and hole cylinders. Here we show that a SF 
pairing interaction favoring exactly such a state is present 
in this material and we discuss the expected consequent 
physical properties. 

The SF spectrum is unusually rich for this compound 
and comes from three separate sources. First, the sys- 
tem is relatively close to a Stoner ferromagnetic instabil- 
ity. Second, there is a nearest-neighbor antiferromagnetic 
(AFM) superexchange [2fl| (d) . Third, there are nesting- 
related AFM spin-density-wave type SF near wave vec- 
tors connecting the electron and hole pockets. The latter 
appear to be the strongest ones. The corresponding in- 
teraction connects the well-separated FS pockets located 
around T, and around M. Though repulsive in the sin- 
glet channel, these would nevertheless be strongly pairing 
provided that the OPs on the two sets of the FSs have 
opposite signs. The main message of our paper is that 
this "s± " superconducting state is both consistent with ex- 
perimental observations and most favored by SF in this 
system. 



As opposed to the undoped material 3, we do not 
find any FM solution for the doped compound, even 
with the more magnetic GGA functional. This sug- 
gests that the main function of doping is to move the 
system away from a ferromagnetic instability (see, how- 
ever, Ref. [20|)- That the system becomes less mag- 
netic is due to the fact that upon doping the heavy 3D 
hole pocket near T rapidly fills and the total DOS drops 
(by a factor of two). At a doping level x = 0.1, we ob- 
tain in the GGA an unrenormalized Pauli susceptibility 
Rexo(q = 0, lj = 0) w 4 x 1(T 5 emu/mole (A(0) = 0.64 
states/eV/spin/Fe). Using 1.1 eV for the Stoner I on Fe, 
we obtain a renormalized %(0) of 0.14xl0 -3 emu/mole, 
a large renormalization, but much smaller than what is 
needed to explain the experimental value [lj]. Note that 
in the undoped system the calculated susceptibility is 
larger, and the experimental one smaller, than in the 
doped one. This suggests that besides FM SF there are 
other, more important spin excitations in the system. 

The first candidate for these is a superexchange 
corresponding to the standard simple nearest-neighbor 
checkerboard antiferromagnetism. Importantly, there is 
also substantial direct Fe-Fe overlap |8j, which leads to 
an additional AFM exchange of comparable strength and 
with the same checkerboard geometry. 

Further discussion requires an understanding of the 
fcrmiology in clearer terms. The band structure may 
at first appear intractably complex, but it is in fact rel- 
atively simple, involving only three Fe orbitals near the 
Fermi level. The hole pockets around T originate from Fe 
d xz and d yz bands that are degenerate at T, and form two 
nearly perfect concentric cylinders. The electron surfaces 
are better understood if we recall that the underlying Fe 
layer forms a square lattice with the period a = a/y2 
(Fig. [2|). The corresponding 2D Brillouin zone (BZ) is 
twice larger and rotated by 45°. If we could "unfold" the 
Fermi surface of Fig. [TJ we would find the same two hole 
pockets at T and one electron pocket around the point X 
in the large BZ. The latter is formed by the d yz band (or 
d xz near Y) that starts 180 meV below the Fermi energy, 
disperses up along X M and is practically flat along X 
r. This band is hybridized with the d xy band (or d x 2_ y 2 
in the two-Fe cell), which starts from X at an energy of 
-520 meV, and is instead dispersive along X T . Upon 
hybridization, these two bands yield an oval cylindrical 
electron pocket, elongated along TX [111 ]. 

Electronic transitions from the hole pockets to the elec- 
tron pocket should lead to a broad (because the electron 
pocket is oval rather than circular) peak in the nonin- 
teracting susceptibility xoiQi 1 - ~ * 0), at q = (7r/5,0), 
while the superexchange interaction J(q) on the square 
lattice should be peaked at q = (tt/o,, tv/cl) (T in the 
downfolded BZ). The renormalized susceptibility, x(q) = 
Xo (?)/[! — J(o)Xo(q)] then has a rich structure with max- 
ima at F, X and M. For the true unit cell with two Fe, 
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FIG. 3: color online The imaginary (a) and the real 
(b) parts of the non-interacting susceptibility xo{q,ui — » 
0), in arbitrary units. Within common approximations 
Imx=Im[ ^ ^ ] is measurable by neutron scattering 

and Rex controls the pairing interaction, in the sin glet chan- 
nel proportional to 1/[1 — J(q)xo(q,0)] (see Ref. [13fl for a 
review.) Note that the RPA enhancement will strengthen 
both peaks. 



both f and M fold down into the T point, while X folds 
down into the M point. The corresponding folding of the 
Fermi surfaces makes the electron pockets overlap, form- 
ing two intersecting surfaces (Fig. [T]). It is important 
to appreciate from this gedanken unfolding that already 
on the level of the noninteracting susceptibility there is a 
tendency for antiferromagnetic correlations with a wave 
vector different from the superexchange one. 

In Fig. [3] we present the calculated [T3| xo(<3S w ) = 
f^-e^-w^iS a t w ~ * and q z — n/c (x is practically 
independent of q z ), The peak at M, derived from inter- 
band transitions, is very broad, as expected, with some 
structure around the M point, due to the particular orien- 
tation of the two oval pockets at M and the size difference 
between the hole and electron cylinders for finite doping. 
With minor modification, this structure is present also in 
the undoped compound. 

We have also performed magnetic calculations (dis- 
cussed in a separate publication) in a supercell corre- 
sponding to the superexchange, q = M (=T in the down- 
folded Brillouin zone) and nesting-induced q = X (=M) 
spin density waves. We find, indeed, that the tendency to 
ferromagnetic ordering is suppressed in the doped com- 
pound, while the tendency to nesting-based antiferro- 
magnetism is present, even leading to an actual insta- 
bility at the mean field level. It is worth noting that 
while this instability does not appear in actual supercon- 
ducting materials, it has been now found experimentally 
in the low doping regime [2fj| . 

The strong AFM SF around M favor our proposed s± 
state. Cases where SF-induced interactions connect two 
pockets of the Fermi surfaces, including SF originating 
from electronic transitions between the very same pockets 
have been considered in the past. 14] However, they 
involved FS pockets related by symmetry, which strongly 
restricts the phase relations between them. In our case 
the two sets of pockets are not symmetry related, and 



nothing prevents them from assuming arbitrary phases. 

For the singlet case, the coupling matrix between the 
hole (h) and the electron pockets (e) is negative, X e h < 
0. The diagonal components emerge from competition 
between the attractive phonon-mediated and repulsive 
SF-mediated interactions and are, presumably, weak. If 
A?t is the average of the phonon-mediated interaction 
over the wave vectors q < < 2k F , and X s J h the same for 
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X S ff h . Our calculated electron- 



the SF, then Xhh = X 
phonon coupling comes mostly from small wave vectors, 
that is, \l h h +\*£ ~ 0.2, and A^ 1 » 0. Phonons, therefore, 
though weak, promote the s± state. On the other hand, 
X S f/ h and Ag{, come from FM (small-q) SF and are pair 
breaking. Finally, the superexchange interaction does 
not affect Xhh, but the corresponding wave vector q = 
TM connects the e pockets near X and Y and are also 
pair breaking 20] (c). 

The T c , as usual for a two-gap superconductor [l5j] . 
is defined by the maximal eigenvalue of the A matrix, 
•"W/ = [Xhh + A ee + ^4A e / i A/ ie + (A ee - Xhh) 2 ]/ 2, and 
the OPs A e ,/t are defined by the corresponding eigen- 
vector. In our case, the signs of A/, and A e will be 
opposite, and their absolute values (despite the differ- 
ent densities of states) will be similar since presumably 
(A ee — X h h) 2 <C X e hXhe- This "s±" state is analogous to 
the states proposed previously for semimetals(l6j and bi- 
layer cuprates[17|. 

From the point of view of neutron scattering, the struc- 
ture of the peak in x near M is important, but for SF 
induced superconductivity it does not matter at all. A 
well defined spin-excitation requires a sharp peak, but the 
pairing interaction is integrated over all possible q vec- 
tors spanning the two sets of Fermi surfaces so that only 
the total weight of the peak is important. Dimensional- 
ity is however very important, as 3D coupling supports 
a long range magnetic ordering (competing with super- 
conductivity). Indeed, the experiment shows an abrupt 
appearence of superconductivity at x « 0.03, roughly 
where the 3D pockets around T disappear. 

One might envision a triplet state similar to the de- 
scribed singlet one, fully gapped, as expected for uni- 
tary 2D p-wave states, and with different amplitudes (or 
signs) on the two cylinders. The similarity, however, is 
misleading. In the triplet channel, SF induce attrac- 
tion, but given the relatively large width of the AF peak 
(FigEl), a large part of the pairing will be lost as only 
SF with a wave vector exactly equal to (n/a, iv/a) will be 
fully pairing, and some others will even be pair-breaking. 
Therefore, in this scenario where the antiferromagnetic 
SF around the M point provide the primary pairing in- 
teraction, we expect the lowest energy superconducting 
state to be s±. The structure of the OP in real space can 
be evaluated using the lowest Fourier component, com- 
patible with the proposed s± state, namely cos k x +cos k y 
(equivalently, cosfc^ cosfcy). This corresponds to pairing 
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of electrons that reside on the nearest neighbors, just as 
in the d-wave case, so that the onsite Coulomb repulsion 
is not particularly destructive for this state. 

Finally, we discuss the experimental ramifications of 
the s± state. In many aspects the y a re similar to that 
of the s± state proposed in Ref. [l7| and discussed in 
some detail in Refs.[l7, 18, The thermodynamic and 
tunneling characteristics are the same as for a conven- 
tional two-gap superconductor The two-gap charac- 
ter, however, may be difficult to resolve, given the dom- 
inance of the interband interactions which will render 
the two gap magnitudes similar. Nonmagnetic small-q 
intraband (e — e or h — h) scattering, as well as the in- 
terband spin-flip pairing will not be pair breaking, but 
paramagnetic interband scattering will, resulting in fi- 
nite DOS below the gapQ, consistent with specific heat 
measurements [5| . The most interesting feature of the s± 
state [n| is that the coherence factors for exciting Bo- 
golyubov quasiparticles on FS sheets with opposite signs 
of the OP are reversed compared to conventional coher- 
ence factors. We outline a few important consequences 
of the relatively straigtforward application of this con- 
cept to experimental probes. First, one expects a qual- 
itative difference between experiments that probe verti- 
cal transitions (q=0) and those that probe transitions 
with q close to 7r/a,7r/a. For instance, the spin suscep- 
tibility at q=0 will behave conventionally, i.e. exponen- 
tially decay below T c without any coherence peak, while 
the susceptibility for qw ir/a,ir/a will have a coherence 
peak that should be detectable by neutron scattering as 
an enhancement below T c jl9| . AFM SF near 7r/a,7r/a 
dominate in the doped material, and so the usual coher- 
ence peak in the NMR relaxation rate, which averages 
equally over all wave vectors is expected to d isap pear 
or be strongly reduced. It was shown in Ref. [l8j that 
Josephson currents from FSs with different signs of the 
OP interfere destructively, and the net phase corresponds 
to the sign of the FS with the higher normal conductance. 
In the constant relaxation time approximation, both in- 
plane and out of plane conductivities are dominated by 
the electron pockets. This is unfortunate, since there 
would otherwise be a tt phase shift between the ab and c 
tunneling and corner-junction experiments could be used, 
as in high-T c cuprates. 

To summarize, we argue that the fermiology found in 
doped LaAsFeO gives rise to strong, but broad antiferro- 
magnetic spin fluctuations near the M point in the Bril- 
louin zone, while the tendency to magnetism existing at 
zero doping is suppressed. These fluctuations, while too 
broad to induce a magnetic instability, are instrumental 
in creating a superconducting state with OPs of opposite 
signs on the electron and hole pockets. 
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